function [Matrix,RHS] = DLS_GSBuilder(M,Order)
%Build matrix and RHS
deg=Order+1;
n=length(M);
Matrix=zeros(n,deg);
RHS = M(:,2);
for k=0:deg-1
Matrix(:,k+1)=M(:,1).^k;
end
end